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Three-loop QCD corrections and 6-quark decays 

John Paul ArchambaullQ and Andrzej Czarnecki 
Department of Physics, University of Alberta 
Edmonton, AB, Canada T6G 2J1 

We present three-loop (NNNLO) corrections to the heavy-to-heavy quark transitions in the limit 
of equal initial and final quark masses. In analogy with the previously found NNLO corrections, 
the bulk of the result is due to the /3o a s corrections. The remaining genuine three- loop effects for 
the semileptonic b —* c decays are estimated to increase the decay amplitude by 0.2(2)%. The 
perturbative series for the heavy-heavy axial current converges very well. 
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The b quark decays most frequently into another relatively heavy quark, the c. These decays, especially when 
semileptonic, provide insight into the properties of heavy hadrons and an opportunity to determine Standard Model 
parameters such as the Cabibbo-Kobayashi-Maskawa (CKM) matrix element V c b 0. 

In order to achieve good accuracy in those studies, one needs to know to what extent gluon emissions modify the 
weak transition amplitude. The evaluation of such effects is complicated by the presence of three energy scales: the 
masses of the radiating b and c quarks, and the momentum transfer to the leptons emitted in the decay process. 
The effects of a single gluon (next-to-leading order, NLO) can be determined as an analytical function of all three 
scales (see e.g. an d the references cited there), but in the next-to- next-to-leading order (NNLO), an analytic 

expression is known only in the limiting case in which the c quark remains at rest with respect to the b quark 0, 
f^S ' known as the zero-recoil configuration. It corresponds to the emission of the charged lepton and the neutrino in 
"p^ | opposite directions. 

Beyond the zero-recoil limit, numerous expansions have been employed to study semileptonic b decays with NNLO 
accuracy 0,0,13- Those results improved our knowledge of the b — > c decay rate, the CKM parameter V c b, and in 
Oh. the future may also contribute to the determination of V u b- 
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FIG. 1: Examples of three-loop corrections to the semileptonic b — » c decay. The W can also be emitted from the quark line 
in the points denoted by crossed circles. This figure shows only examples of abelian corrections, similar to photonic vertex 
corrections. The complete list of such diagrams is given in Q. Non-abelian diagrams are shown in subsequent figures. 



In the present paper we study effects of the next-to-next-to-next-to-leading order (NNNLO) to the transition of one 
heavy quark into another. Effects in this order, ©(af), have never before been calculated for any charged particle 
decay. The main reason is that the NNNLO calculations, especially with massive particles in intermediate states, are 
very difficult. 

Here, we consider a particular kinematical limit in which the masses of the initial and final state quarks are equal. 
This is the extreme zero-recoil limit: there is no momentum transfer into leptons and no real gluon radiation. The 
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zero-recoil configuration can be realized also with unequal quark masses, when the leptons, emitted anti-parallel to 
each other, carry away all released energy. Radiative corrections in such configuration also consist of virtual gluon 
exchanges only. If the mass difference between the quarks is not very large, loop diagrams can be evaluated by 
expanding in this mass difference, around the equal mass ("extreme zero-recoil") limit considered in this paper. 

Taking equal quark masses greatly simplifies the NNNLO calculation because we deal with two-point Green func- 
tions, similar to a quark propagator, with zero-momentum insertions (see Fig. [Q. In fact, the technically related 
problems of the three-loop anomalous mag netic moment, Dirac form-factor, and heavy quark renormalization con- 
stants have already been solved 0, 0, Il2j | . 

For the present project, we modified a program originally written for the latter two studies [l3| . Thus, with 
relatively little effort, we can obtain the first NNNLO correction to a heavy-to- heavy quark transition and confirm 
the expectations of its moderate size. It should be mentioned that despite the special kinematic limit in which we 
work, our result may help estimate NNNLO effects in a variety of other configurations because most of the kinematic 
effects reside in easily determined phase space effects, while the matrix elements depend on kinematics only relatively 
weakly. (There are exceptions to this picture; for example, along the zero-recoil line, the matrix elements are affected 
by a mass-singularity effect, manifesting itself in logarithms of the initial- and final-state quark mass ratio |l4l Il5j|. 
However, the physical ratio of b and c quark masses is not very large and our result obtained in the equal mass limit 
is expected to hold with reasonable accuracy |l6(.) 



II. CORRECTIONS TO THE HEAVY-TO-HEAVY TRANSITION b -» civ 



The effect of QCD corrections on the Web vertex is parameterized using two functions, 

7m(! ~ 7s) -> 7m [vv(q 2 ) - VA{q 2 )ls\ ■ (1) 

In the extreme zero-recoil limit we are considering, the momentum transfer into the W vanishes. In this case, the 
vector part of the interaction does not receive any QCD corrections, 

r, v (q 2 = 0) = 1. (2) 

This is a useful check: the vector coupling part of the sum of all our three-loop diagrams, after the quark wavefunction 
renormalization, must vanish. 

The axial part of the interaction does receive finite corrections even in the q 2 — limit. They are expressed as a 
power series in a Sl 

VA (q 2 = 0) = ru = 1 + ^C F ^ + C,r,$> + C FV f + Otf) . (3) 

Throughout this paper, we define a s in the MS scheme at the scale of the quark mass, and use the pole definition of 
the quark mass. 

The first two corrections to t)a have been known for a long time 0, IT7L Il8| , 
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V A - 2 , 

(2 ) ( 373 tt 2 \ , ( 143 7T 2 7r 2 ln2 1 \ /115 1 ,\ 7 

174 = C F [ + — + (C A - 2C F ) + Cs + N H T R 7T 2 + — Nl.Tr. 

IA *l 144 6 / 1 7 V 144 12 6 4 Sc 7 U 36 3 / 36 



(4) 



The SU(3) factors are C F — 4/3, Ca = 3, and Tr = 1/2, and Nh,l denote the numbers of heavy and light quark 
flavors. For b — > c decays, Nh — 2 and Nl — 3. The two- and three-loop corrections are expressed in terms of the 
Riemann zeta function and a tetra-logarithm, 

C 3 ^ 1.202056903, 
C 5 ^ 1.036927755, 

a 4 = Li 4 Q^j ~ 0.5174790617. (5) 
The result of the study reported here is the third correction, 77^ . We find 

ia I 162 27 
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Among these ten contributions, differing by SU(3) factors, the first seven come from diagrams containing one or two 
fermion loops, with light and/or heavy quarks. Among them, the second and the fifth arise from non-abelian diagrams 
with three- or four-gluon vertices, shown in Fig. [21 The remaining five contributions containing fermion loops are 




FIG. 2: Three-loop non-abelian diagrams with fermion loops. Solid line closed loops denote light or heavy quarks. As in the 
previous figure, the W boson can be emitted at any vertex denoted by a crossed circle. 



analogous to those in QED, with one example shown in Fig.Q] 

The last three contributions in Eq. (j^J correspond to diagrams without closed fermion loops. The first and the 
second of them are non-abelian, and receive contributions from the diagrams shown in Fig. as well as from the 
non-planar diagrams without multi-gluon vertices. The latter, QED-like diagrams, also give the last contribution in 
Eq. ©. 

In addition to the vertex correction, the wavefunction renormahzation also contributes to 77^. The O(af) part of 
the on-shell quark wavefunction renormahzation constant Z2 was found in To our knowledge, the present study 
is the first application of that result. 

One novel feature of the on-shell constant Z <? at the NNNLO is that some of the non-abelian contributions are 
gauge-dependent in dimensional regularization [l2|. In contrast, dimensionally regularized on-shell Z2 in QED is 
gauge-independent to all orders, and before the explicit calculation in 12], it had been conjectured that the QCD 
result might also be gauge-independent 0]. The calculation presented here shows that the gauge dependence of 

Z2 found in |T^ is needed to make the full vertex correction tyf v gauge invariant. The cancellation of the gauge 

(3) 

parameter dependence is a useful check, in addition to the vanishing of rjy. 
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FIG. 3: Three-loop non-abelian gluonic diagrams. 



III. NUMERICAL RESULTS AND SUMMARY 

Substituting the numerical values of the coefficients of the SU(3) factors, we find for the first three corrections to 

rff = -0.5, 

r^ 2) ~ -0.976 Ca + 1.01 C F - 0.0954^^ + 0.194^^, 

r^ 3) ~ -2.06 C\ + 0.958 C A C f - 2.15 C% - 0.124 C A N H T R + 1.03 C F N H T R + 2.13 C A N L T R 

+0.202 C F N l T r - 0.0507 N% T 2 R - 0.0647 N H N L T R - 0.288 N 2 L T%. (7) 

For the case of the b quark decay into a c quark (Nl = 3, Nh = 2), we get in the extreme zero-recoil limit, 

VA ~ 1 - 0.667^- 1.85 (— ) 2 - H.l (— ) 3 + 0(«l) 
~ 1 - 0.0510- 0.0108- 0.00495 

~ 0.933. (8) 

In the last two lines we substituted a s — 0.24, roughly (to within 20%) corresponding to the value appropriate for 
the scale ^mt,m c . We see that the three-loop calculation presented here contributes about 7% of the deviation of t]a 
from unity, and decreases t]a by about half of a percent. A large part of this NNNLO correction is due to the known 
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effects that can be absorbed in the running of the coupling constant, that is contributions of the type oP 8 0q, with 
(3 = 11 — ^Np, and Np denoting the number of fermions contributing to the running. This is often referred to as the 
Brodsky-Lepage-Mackenzie (BLM) corrections [2(J. We can find it using the coefficient of (diagrams containing 
two light fermion loops) from Eq. JSJ. We obtain, taking Np = 4 quarks in the running, 

^<?><BLM) = ( 4 -f) 2 C^(JL-!!) 

= -15.0. (9) 

We see that the remaining, "genuine" , three-loop correction is small, 

C^ 3) (non-BLM) = -11.1 + 15.0 = 3.9. (10) 

It is possible to extend this result beyond the limit of zero momentum transfer to the leptons by expanding in the 
difference of the heavy quark masses. However, we believe that our result estimates the true value of the correction 
at physical values of uib and m c to within at most a factor of two. The corrections due to unequal m,b and m c are 
suppressed by two powers of 1 — m c /mi, 16]. Given that the new non-BLM effect in Eq. I|10|) increases t]a by about 
0.2%, we believe that it is safe to estimate the corresponding effect for the real b — > c decay at zero recoil as 0.2(2)%. 
Finally, we note that t\a determines corrections to the total semileptonic decay rate if the difference of the quark 
masses is relatively small [2lJ and can be used to estimate the b — > c decay rate 22] . Such analysis will be presented 
elsewhere. 
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